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Abstract 

SU(N) Lie algebras possess discrete symmetries which can lead naturally to 
stable vector dark matter (DM). In this work, we consider the possibility that 
the dark SU(N) sector couples to the visible sector through the Higgs portal. We 
find that minimal CP-conserving hidden “Higgs sectors” entail stable massive 
gauge fields which fall into the WIMP category of dark matter candidates. For 
SU(N), N > 2, DM consists of three components, two of which are degenerate 
in mass. In all of the cases, there are substantial regions of parameter space 
where the direct and indirect detection as well as relic abundance constraints 
are satisfied. 
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1 Introduction 

The puzzle of dark matter remains an outstanding problem of particle physics. One of 
the more attractive approaches to this problem exploits the fact that weakly interact¬ 
ing massive particles (WIMPs) in thermal equilibrium produce the relic dark matter 
abundance in the right ballpark. Particles of this type appear in many extensions of 
the Standard Model. In this work, we explore the possibility that the Standard Model 
is connected to the dark sector through the Higgs portal [T], 

4"portal = , (1) 

where H is the Higgs field and 0 is the “hidden Higgs”, that is, the field responsible 
for breaking the gauge symmetry of the hidden sector. We assume the dark sector 
to be endowed with U(l) or SU(N) gauge symmetry. In that case, the massive gauge 
fields can play the role of WIMP-type dark matter quite naturally. Indeed, they are 
weakly coupled to the Standard Model and, owing to inherent discrete symmetries, 
can be stable. 

The U(l) case was considered in Ref. [2] where the stabilising Z 2 symmetry was 
found to be related to charge conjugation. The SU(2) example was worked out in 
Ref. [3] and the DM stability was attributed to a custodial SO(3). In this work, we 
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extend these ideas to larger SU(N) Lie groups and uncover the nature of the underlying 
discrete symmetries. We focus on weakly coupled theories, although conhning hidden 
sectors represent a viable alternative HE]. Finally, we perform a comprehensive 
study of DM phenomenology in all of these cases, which includes direct and indirect 
DM detection as well as an analysis of the DM relic abundance. Our main conclusion 
is that massive hidden gauge fields serve as viable and attractive DM candidates. 

2 Massive U(l) and SU(2) gauge fields as vector 
DM 

In this section, we review the cases of massive U(l) and SU(2) gauge fields as vector 
DM and identify the underlying symmetries leading to stability of vector dark matter. 
In what follows we assume that the hidden sector consists only of gauge fields and of 
scalar multiplets which are necessary to make these gauge fields massive. 

2.1 Hidden U(l) 

An Abelian gauge sector provides the simplest example of the vector DM model en¬ 
dowed with a natural Z 2 symmetry j2]. In this case, the Z 2 corresponds to the charge 
conjugation symmetry. 

Consider a U(l) gauge theory with a single charged scalar 0, 

-Chidden = > (2) 

where is the held strength tensor of the gauge held and D(0) is the scalar 
potential. For easier comparison to the non-Abelian case, we take the charge of 0 
to be -1-1/2. Suppose at the minimum of the scalar potential 0 develops a VEV, 
(0) = l/\/2 V . The imaginary part of 0 gets eaten by the gauge held which now 
acquires the mass ttia = gv/2, where g is the gauge coupling. The real part of 0 
remains as a degree of freedom. Denoting it by p and normalising it canonically, 
0 = 1/v^ (p + h) , we get the following gauge-scalar interactions: 

A^A^ . (3) 

The system possesses the Z 2 symmetry 

) ( 4 ) 

which is the usual charge conjugation symmetry. In terms of the original scalar held, 
this symmetry acts as 0 —)■ 0* and —>■ — A^, which is preserved by both the 

Lagrangian and the vacuum. The Z 2 makes the massive gauge held stable. Note that 
this symmetry applies to a sequestered U(l) which has no tree level mixing with the 
hypercharge, in which case no mixing is generated radiatively either. 

Interactions with the visible sector proceed through the Higgs portal coupling 

-^^portal = —^h(l>\H\^\(j)\^ , (5) 
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which also leads to the Higgs mixing with p. In the unitary gauge, the Higgs held is 
given by = {0,v + h)/ a/ 2. The helds p and h are then to be expressed in terms of 
the mass eigenstates hi ,2 as follows: 

p = —hi sin 6 + h 2 cos 6 , 

h = hi cos 9 + h2 sin 6* , (6) 

where the mixing angle 6 is constrained by various experiments, most notably, by LEP 
and LHC. The upper bound on sinh depends on the mass of the heavier state ^2 and 
is around 0.3 for of the order of a TeV, see e.g. hgure 3 of Ref. [6] for details. The 
lighter state hi is identihed with the 125 GeV Higgs, while the mass of the second 
state can vary in a wide range. 

As the Higgs portal necessarily preserves the Z 2 symmetry, is a viable DM 
candidate. All the relevant scattering processes proceed through an exchange of hi 
and h 2 , which include DM annihilation into the SM particles as well as DM scattering 
on nucleons. 


2.2 Hidden SU(2) 

The U(l) considerations can easily be extended to SU(2), albeit with a modihcation 
of the stabilising symmetry. Consider an SU(2) gauge theory with one doublet 0, 

^hidden = >“■ + (D„4,)'D^4, - V(4,) , (7) 

where a = 1,2, 3. The potential is assumed to have a minimum at a nonzero VEV of 
(f). In the unitary gauge, / takes the form 


0 = 


1 



( 8 ) 


with p being a real field and v being the VEV. Denoting the gauge coupling by this 
leads to the gauge boson mass tua = 9^/2. The scalar-gauge and gauge-gauge held 
interactions are given by 

a;a‘ “ + “. 

ADg.g = <7)2 _ _ ^g) 

Although the triple gauge vertex breaks the parity of the previous section, it follows 
that the system possesses a. Z 2 y- Z'^ symmetry, 

Z2 : -)■ -A^ , A^ -)■ -A^ , 

^'2'- ■ ( 10 ) 


As a result, all three A“ helds are stable and can play the role of dark matter. While 
the above symmetry is sufficient to ensure stability of DM, it actually generalises in 
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this simple case to a custodial SO(3) |3]. As we will see, for larger SU(N) groups, 
it is the discrete symmetry that plays a crucial role. The hrst Z 2 is associated with 
a gauge transformation, while the Z '2 generalises the charge conjugation symmetry, 
i.e. it corresponds to complex conjugation of the group elements. 

As before, the dark sector couples to the visible one through the Higgs portal. 
The discussion of the Higgs mixing with p of the previous section applies 
here as well. Clearly, the Z 2 x Z'^ symmetry is preserved by the Higgs portal and the 
hidden gauge helds couple to the visible sector only in pairs. 

3 Extension to SU(3) 

While the SU(2) case is straightforward, larger SU(N) groups exhibit more compli¬ 
cated breaking patterns. In a phenomenologically viable set-up, the symmetry must 
be broken completely to avoid the existence of massless fields (barring confinement). 
One possibility would be to break SU(3) by two scalar multiplets in the fundamen¬ 
tal representation, i.e. triplets. One may also explore other options involving SU(3) 
tensors. As we detail below, our conclusion is that a single irreducible representation 
with two indices cannot break SU(3) completely, which leaves the two triplet option 
as the minimal one. 

3.1 Breaking SU(3) by tensor fields 

The lowest order SU(3) tensor whose generic VEV can break SU(3) completely is the 
symmetric tensor that is 6 of SU(3). Gauge transformations act on (pij as 

(f)^U(t)U^ . ( 11 ) 

By virtue of Takagi’s matrix decomposition, this allows one to bring (pij to the diagonal 
form, 

/0i 0 0 \ 

0=0 02 0 , ( 12 ) 

yo 0 03 / 

where 0i, 02,03 are real up to an overall complex phase. If the VEVs of 0i, 02 and 03 
are all different, SU(3) is broken completely. However, when some of them coincide, 
the residual gauge group is at least SO(2). 

In order to determine what VEVs are possible, let us write down the most general 
gauge invariant potential for 0jj. It is easy to convince oneself that the potential has 
the form 


V = m^ Tr0V + AiTr(0V)^ + A 2 (Tr0V)^ + (/W Det0 + h.c.) , (13) 

where m? can be negative. The minima of this potential determine the 0, VEVs. 
Acting upon V with the operator 0j0/00j (no summation over i), one hnds that all 
nonzero (00 satisfy the same equation. This implies that (00 are either degenerate 
or zero. For the case p = 0, such behaviour has been noticed in [TJ IB]- We find 
(analytically and numerically) that it persists in the case of nonzero p as well. As 
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a result, the residual gauge symmetry is at least S0(2) and the corresponding gauge 
bosons remain massless]^ Thus, the model with a single symmetric tensor is unrealistic. 

One may also consider the possibility of breaking SU(3) by an antisymmetric ten¬ 
sor. By virtue of Youla’s decomposition, it can be gauge-transformed to the block- 
diagonal form 

/ 0 01 0 \ 

0 = -01 0 0 . (14) 

\ 0 0 0 / 

Therefore, the residual gauge symmetry is at least Sp(l)=SU(2), which can also be 
understood via the equivalence between the antisymmetric tensor and the (anti)funda- 
mental representation. Again we find that the model is unrealistic. 

Similarly, complete SU(3) breaking cannot be achieved by a VEV of an adjoint 
scalar. In this case, the group rank is preserved and thus there are massless gauge 
helds. 

The next simplest option is to combine the symmetric and antisymmetric ten¬ 
sors. This system would have 9 complex degrees of freedom, while two SU(3) triplets 
only have 6 . Thus SU(3) breaking with two triplets is minimal and sufficient for our 
purposes. 


3.2 Breaking SU(3) by triplets and x 


Misaligned VEVs of two triplets break SU(3) completely. This breaking pattern can 
be understood in stages: the first triplet VEV reduces the symmetry to SU(2), while 
the second breaks the remaining SU(2). VEVs misaligned in SU(3) space represent 
a generic situation, that is, they result from the minimisation of a general scalar 
potential consistent with SU(3) symmetry. Therefore, such a breaking pattern is 
phenomenologically viable. 

Before going into details, let us identify the Lie group discrete symmetries which 
eventually lead to DM stability. One way to find them is to analyse the SU(3) structure 
constants (using the usual Gell-Mann basis). 


/ 

/ 

/ 


123 

147 

458 


1 , 


_^156 
^678 ^ 



= f 


345 


-/ 


367 


1 

2 ’ 


( 15 ) 


Identifying the transformation properties of the generators with those of the gauge 
fields, we define the “parity” transformation as 




It is easy to see that the structure constants are invariant if the parities are 


( 16 ) 


Z 2 '■ h(®) = —1 for 0 = 1) 2,4, 5, 

T]{a) =+1 for a = 3, 6 , 7, 8 , (17) 

^ Ref. [1] has considered (0) proportional to the unit matrix, which entails unbroken SO(3). The 
corresponding gauge bosons may be confined in glueballs. 
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and also 

Z 2 : r]{a) =—1 for a = 1,3,4, 6 , 8 , 

r]{a) =+1 for a = 2, 5, 7. (18) 

One may notice that the hrst Z 2 acts on the off-diagonal generators with nonzero 
entries in the hrst row, while the second rehects the real SU(3) generators. In Section]^ 
we will show that these symmetries generalise to arbitrary SU(N) and that the hrst Z 2 
is a gange transformation, whereas the second corresponds to an onter antomorphism 
of the gronp, i.e. complex conjngation of the gronp elements. 

These symmetries are inherited by the Yang-Mills Lagrangian. If CP is conserved, 
they are also preserved in the matter sector leading to stable dark matter. Below we 
study the relevant interactions in detail. 


3.3 Explicit example 


Consider an SM extension by two complex helds and 02 transforming as triplets of 
hidden gauge SU(3). The Lagrangian of the model is 


C 


SM 


c 


portal 


+ c 


hidden 


( 19 ) 


where 

-CsM D CsM = ^\H\^ + ml\H\^ , (20a) 

—-Cportal = hportal = + ^H22 |77|^|02|^ — {^H12 \H\'^4’l4>2 + h.C.) , (20b) 

'^hidden = ~ ~ 14idden • (20c) 

Here, = d^Ay — dyA^ + ig[A^y^ Ay] is the held strength tensor of the SU(3) gauge 

helds D^(j)i = d^(l)i + igA^(l)i is the covariant derivative of 0*, H is the Higgs 

doublet, and the most general renormalisable hidden sector scalar potential is given 
by 


Vhidden(01, 02) = "in|0lP + "^22102^ “ (mi20l02 + h.C.) 

+ + “^| 02 |"^ + -^ 31011 ^ 1021 ^ + -^ 4101021 ^ 


2 ' 2 

t^_^2 , \ AAA’i/A \ J2/'M 


-^5 


“^(0102) +A6|0i| (0102) + -^7|02| (0102) + h.C. 


( 21 ) 


Using SU(3) gauge freedom, 5 real degrees of freedom of 0i and 3 real degrees of 
freedom of 02 can be removed. Therefore, in the unitary gauge 0i, 02 read 


01 = 


1 

71 


0 

0 


Vl + (fi 



V2 + (^2 

(^3 + Pz) +i{VA + </54) 


( 22 ) 


where the Vi are real VEVs and (pi-^ are real scalar helds. Analogously, we express 
the Higgs held in the unitary gauge as = (0,n + h)/y/2. 

In what follows, we make two assumptions which are crucial for stability of vector 
dark matter: 


6 




• the scalar potential is CP invariant 

• the VEVs of 01 , 02 are real. 

The hrst condition implies that the scalar couplings are real, while the second assumes 
that no spontaneous CP violation occurs (^4 = 0). As a result, C'P-even and C'P-odd 
helds do not mix. 

In this case, the Z 2 x Z 2 symmetry extends to the Higgs sector as well. Under 
the hrst Z 2 all are even, while the second Z 2 rehects <^94 and leaves the other helds 
intact. As we detail in Section]^ this assignment follows from the explicit form of the 
hrst Z 2 as a gauge transformation and the fact that the second Z 2 acts as complex 
conjugation. In any case, these are explicit symmetries of the Lagrangian and the 
vacuum. The full list of the parities is presented in Table Clearly, the lightest 
states with non-trivial parities cannot decay to the Standard Model particles. 


helds 

Z2 X Z2 

h, (p0 

(+) +) 

42 45 

(“! +) 



<^^A^,A®,A® 

(+) “) 


Table 1: Z 2 Z 2 parities of the scalars and dark gauge bosons. 


We now discuss the Lagrangian in more detail, starting with the covariant derivates 
of 01 and 02 , 

|D^0,|2 = |a^0,|2 + i~gAl ((a^00tT“0, - h.c.) + . (23) 


Inserting the parametrization (22) with V 4 = 0, we get the kinetic terms for the scalars 
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(24) 


Z=1 


the mass terms for the gauge helds, the mixing terms as well as the gauge-scalar 
interactions. Let us hrst discuss the terms quadratic in the helds. The third term on 


the r.h.s. of Eq. (23) contains the mass terms for the gauge bosons. 


^ ^ 2 ("^b4)-^(l,4)Al,4) + ^,5)-^(2,5)A(2,5) + A^,6,8)-^ (3,6,8)A(3,6,8) + , 

(25) 

where we have used the shorthand notation = (A|j,A^) and similarly for the 
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other gauge bosons. The mass-squared matrices are 


q 

A^(1_4) = Ai(2.^) = — 


(2,5) 


V2 V2V3 
V2V3 vf + 


V2V3 

2 I „,2 I „,2 


-Ad (3,6,8) = -J I -V 2 V 3 Vf+V^ + vl 


-uI/a/S 

-V 2 V 3 /V 3 


l\pi -V2V3l\pi (4uf + + 4u^)/3 


■^( 7 ) = + vl + vl) . 


(26) 


The second term on the r.h.s. of Eq. (23) contains the A^d^ip mixing terms 


In general, terms of the type KaiA^d^if^' can be removed by the field redefinition 


9 V 3 


9 V 3 


9 V 2 


(27) 


4“ = 4“ + 5„F“ , where F“ = (XF')., km v ' ■ 


(28) 


We then have 


where 


- ttr{G„„G'“'} + i(S„y)'^ + \m^AIA>‘>’ + = 

- yr{G„„G'“'} + I't.MfWV) + \m^a;A'‘'’, 


Gi,u = d^A^ - d^A^ + ig[A^ - d^Y, - d^Y ], 

7q = ^ij ~ Vi^g^M-ab 1 


(29) 

(30a) 

(30b) 


and we have defined Y = Y°‘T°‘. The kinetic terms for the vector fields are still 
canonically normalised, unlike those for the scalar fields. The latter can be normalised 
canonically by a further field redefinition 

(f" = , where (ca^a;)^ = 7 ^-. (31) 

Note that the term — includes additional couplings of the gauge fields to 

the scalars. The resulting vertices are obtained by replacing ^4“ in the triple and quar- 
tic gauge boson terms by —d^Y°‘ such that interactions of the type {A)^dY, {A)‘^{dY)‘^, 
etc. arise. The analysis of the general case is very cumbersome, so further we will 
focus on the simple case of U 3 = 0 which retains all the relevant physics. In that case, 
the couplings involving BY play no role in the DM phenomenology. 


3.4 Detailed study for t ;3 4 = 0 

As long as ui and V 2 are nonzero, SU(3) is broken completely. Hence it suffices to 
consider the case U 3 = 0, which simplifies the analysis. Then the only mixing term 
among the gauge bosons is A^^A^^ and the only gauge-scalar mixing terms are 
Ald^ip3. 

In what follows, we restrict ourselves to the case U 3 = 0. 
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3.4.1 Masses 


Gauge boson masses The gauge boson mass eigenstates are 


4/3 

1 — 


cos a A^, + sin a A^, 




V 

. ' > I where tan 2 a= —^ , 

cos a AZ - sin a AZ } 2 vf - vl 


(32) 


and the masses ar^El 


2 

= —T 


1 - 


tan a 


Vs r 


2 Vvl 1 

m^,8 = 


3 1 


tan a ’ 


(33) 




while the other gauge boson masses can be read off directly from Eq. (26) for ^3 = 0: 


2 2 V 2 2 2 V 2 2 2 V / 2 , 2 \ /oA\ 

tUai = = —^2 , m^4 = m^5 = —v^ , + v^) ■ (34) 


For ^2 < Ui, the light helds are A^ and A'V with the latter being the lightest. It is 
instructive to consider the case v\, so that tana is small and positive. Then A'^ 
is slightly lighter than A^^ and ^4^ by a factor of (1 — while the other five dark 

gauge bosons are all much heavier, by a factor of order vxjv^- The mass degeneracy 
between Aj^ and persists at loop level by symmetry arguments (see Section |^. 

One can easily verify that the heavier states A^"^ and A'^ all decay into the light 
states and the SM particles. The decay proceeds via emission (off-shell or on-shell) 
of the CP even scalars which couple to the SM Higgs and thus to all other SM helds. 

The 3 lightest states all have different Z 2 x Z 2 parities such that they cannot decay 
into each other by emitting SM particles. The only scalar with negative parity is (^ 4 , 
however it is generally heavy (see below) and does not contribute to the above decay. 
Hence, A^ and are stablej^ 


Gauge boson - scalar mixing According to Eq. (28), the A^d^ip mixing terms are 
removed by the redehnition (which does not affect the gauge boson masses) 


Al = Al + , where =- 


V2 4 
- 2 P > 


gvf + vi 


Al = Al + d.Y'^ , where = -- 


V2 


gvf + vi 


2 V, 


(35) 


while, according to Eq. (31), the canonically normalised scalars are 






2Note that ^ + 0 (I), cosa 




= 1 - 


V = 


Vl 


■Yf+ 


-,‘p 


(36) 


3v2 

32vf 


+ O ^(( 1 ^ and sin a = 


VSvj 

Av? 


o 




^The parities allow for a decay of one DM component into two others, however this is forbidden 
kinematically. 
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Scalar masses The general scalar potential has many parameters. To make our 
discussion more transparent, let us assume the symmetry 02 —02 which does not 

affect the essence of our considerations and requires m ^2 = ^hi 2 = Ae = A 7 = 0. In 
this case, it turns out that the potential has no local minima with all of v,vi,V 2 ,vs 
being nonzero such that setting us = 0 is actually required. It should be noted that in 
practice we are considering the limit in which the above quantities are very small but 
nonzero such that the decay channels for heavy particles, e.g. 03, into the SM fields 
are open. 

The VEVs v,vi,V 2 can be expressed in a compact form using the matrix 

( Ah AhiI Ah22 

Ahii Ai A 3 

Ah 22 A 3 A 2 


(37) 


as well as the matrices Aij defined as (—times the matrix obtained by deleting 
the Tth row and j-th column of A, i.e. det Aij is the (i, j)-cofactor of A. One finds 

= —2{rn\j det An + det A21 + m22 det A31)/ det A , 

vl = det A12 + m\i det A22 + ^22 det A32)/ det A , 

v\ = — 2 {mjj det A13 + det A23 + ml2 det A33)/ det A . ( 38 ) 

The mass terms for the scalars are 

— ^ ^ ^ 2 ) ^3 + ^(-^4 — A5)(n^ + vl) 04 , (39) 

where = (h, ipi, 1 P 2 ) and 

/ Xhv"^ XhuVVi Xh22VV2 \ 

= XhuVVi Xivf XiiViV2 . (40) 

\ Xh22VV2 X3V1V2 X2VI ) 

Note that 03 and 04 are generally heavier than since their masses involve 

vi (unless A 4,5 are very small). 

The matrix is positive definite if and only if det A > 0, det A33 = A^Ai — 
A^ii > 0 and Xh > 0. It can be diagonalised by an orthogonal transformation 
O^m^O = diag(m|^^,■' 22 / 73 ), where 


/ C 12 C 13 S 12 C 12 S 13 \ 

O = I —C13C23S12 — S13S23 C12C23 —C 23 'Si 2 'Si 3 + C13S23 1 ( 41 ) 

\ —C23S13 + C13S12S23 —C12S23 C13C23 + S12-S13S23 / 

and we have used the abbreviation Sij = sin= cos9ij. Instead of providing 
the most general formulae, let us focus on a simplified case. Suppose that the ( 12 ) 
and (23) entries of are much smaller than the other matrix elements, in other 
words, that Ann and A 3 are very small. In this case, the Higgs mixing with ip 2 is the 
dominant one. The reason behind this choice is that the DM constituents A^, and 
A'^ all have a significant coupling to (^ 2 , see Table Given that the mixing between 
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the Higgs and ip 2 is substantial, this facilitates DM annihilation into the SM helds. 
Clearly, similar considerations also apply to the general case with all the angles dij 
being significant. 

Assuming small 111 ^ 2 , 11123 , one finds 


6»i2 P 32 mi2 + q , 

023 ~ P21 in23 + 


tan 2013 ~ 


2Xh22VV2 
\ 2 vl - Xhv‘^ ’ 


(42a) 

(42b) 

(42c) 


where Pij = (ni| - niA)/s, q = ni? 3 /s and s = (nif 3)2 + (mfi - ni^ 2 )(m 22 “ ^ 13 ). 
The mass eigenstates are 



Cish — S13932 
Cl3‘^2 + 5130, 


712 


C 13 A 1 

-h 

SlsPl 


723 


5i3</2i 

P 2 

-C 13 P 1 


and the mass-squared eigenvalues are 


~ Ai^i. 


(43) 


(44a) 

(44b) 


The eigenstates hi, 02 are typically the lighter ones, while Ti is heavy (unless Ai is 
very small). In our analysis of DM phenomenology, we retain only the former states. 
In summary, the relevant light fields are the DM components A^, and A'^ as well as 
the mediators hi, ^2 which link the dark sector to the SM fermions and gauge bosons. 


3.4.2 Couplings 

The full list of the couplings is not necessary for our DM studies. The important 
couplings are those with two gauge bosons and one or two scalars at the vertex. In 
terms of the variables A“ and most of these are listed in Table and Table |2 For 
our applications, the couplings of hi and 02 are obtained from these tables using the 


relation (43), in which one may neglect 0 i 2 and 623 - 


We focus on the case Vi S> V 2 so that DM consists of A^, A^ and A'^. Other helds 
with non-trivial parities decay into these states and the SM particles. For instance, 
the processes (^4 —)■ A^A^ + SM and (p^ —>■ A^' + SM are allowed. (When vi and V 2 are 
close, the composition of DM depends on the mass splittings.) DM annihilation and 
scattering proceeds through an exchange of hi and h 2 - Therefore, only the vertices 
involving these helds play a signihcant role. 
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4 Generalisation to arbitrary SU(N) 


SU(N) is broken completely by generic VEVs of iV — 1 fields 0* in the fundamental 
representation. The can be gauge-transformed to the form 


/ 0 \ 


( ° ^ 


( 0 \ 

Pn-1 

0 

, 02 — 

0 

, ••• , 0Af-l — 

0 

Vdi/ 

p? 

(N-2) 

Pn-1 

, (N-l) 

VPiv-1 / 


(45) 


Here the radial helds and the phases are real. We label the scalars such that 
the lightest gauge helds are associated with the SU(2) subgroup which gets broken at 
the last stage by a VEV of We assume that the VEVs as well as the couplings 

in the scalar potential are all real so that CP is preserved in the hidden sector. 

The generalisation of the Z 2 x Z 2 parity to SU(N) is as follows. The transformation 
properties of the gauge helds are identihed with those of the corresponding SU(N) 
generators. The basis of the N{N — 1) oh-diagonal generators can be chosen 

as 


(T°'^)ij — SiaSjb + SibSja , 

{f^% = -^6,JJb + ^S^b6Ja, (46) 

where a = 1,N — 1 and b = 2, ..,N. With the Cartan generators denoted by 
the Z 2 associated with complex conjugation of the group elements acts as 

rjnab _ y _ rjnab r^ab _ y r^ab _ y _ 


This is a well known outer automorphism of SU(N) which entails the corresponding 
symmetry of the Yang-Mills Lagrangian. 

Another Z 2 can be dehned by rehecting the oh-diagonal generators containing 
nonzero elements in the hrst row: 


rjnla _ y _ rj^la 

rpbc _ y rj-ibc 


r^la 

r^bc 


_ y _ r^la 

fbc 


{b,c>2), 


(48) 


It is easy to show that this Z 2 is an inner automorphism. Indeed, it corresponds to 
the group transformation with 

V = diag(—1,1, ...1) . (49) 

The gauge helds associated with the upper left SU(2) block 

ifi=diag(l,-l,0,...,0) (50) 

transform under these parities the same way as did the SU(2) gauge helds under 
Z2 X Z 2 of the previous section. The above gauge transformation of course leaves the 
Yang-Mills Lagrangian invariant. 
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These symmetries are preserved by gauge interactions with scalars in our set-up. 
The Z 2 associated with complex conjugation acts on scalars by reflecting the complex 
phases, which therefore correspond to odd fields under the transformation (47). This 
symmetry is guaranteed by CP invariance of the Lagrangian and preserved by the 
vacuum (assuming no spontaneous CP violation). The second Z 2 is a gauge trans¬ 
formation. On vectors of the form (0, Oi,..., aAr_i), it acts as multiplication by an 
overall constant phase which cancels in all the Lagrangian terms. It is therefore a 
valid symmetry in the broken phase as well. 

As long as the 0* have a zero first component, the interaction vertices contain an 
even number of and T^“. The gauge fields associated with a > 2 are heavier than 
those corresponding to a = 2. By virtue of the vertices involving [k > 2) and 

the matter fields, they decay to the lighter fields such that only the final SU(2) block 
remains stable. (Similar considerations also apply to other heavy gauge fields.) Then 
DM is composed mostly of the aforementioned whose stability is enforced by 
Z 2 X Z'Q 

As in the SU(3) case, DM consists of three components, two of which are degenerate 
in mass, 

m^i/ = m^ 2 / 7 ^ m^3/ , (51) 


where are the mass eigenstates consisting mostly of with some admixture 

of other gauge fields (see the SU(3) example). The degeneracy persists at loop level. 
This can be seen as follows. The SU(N) Lie algebra possesses a discrete symmetry 
which interchanges the real and imaginary generators with nonzero entries in the first 
row. 


rj^la _ y rjn 


la 


rj^la _ y _ rjn 


la 


(52) 


while all the other generators remain intact. This is achieved by the group transfor¬ 
mation 


[/' = e 2?7 diag(f, 1 , ... 1 ) , 


(53) 


which can be recognised as the square root of U in (49). This gauge transformation 
acts on with a zero first entry as an overall constant phase multiplication. Since 
such a phase cancels in all of the Lagrangian terms, (52) remains a valid symmetry 
even in the broken phase. 

Consider now the mass matrix for the gauge fields associated with and T^“. 
By virtue of (48), only fields corresponding to and can mix, while (47) forbids 
a mixing between the tilded and untilded fields. The resulting mass matrix for is 
then identical to that of according to (^. Analogous considerations apply to the 
kinetic terms. Hence the lightest eigenstates have the same mass. 

The resulting DM phenomenology is analogous to that for the SU(3) case. 


"*^As before, we take the “phase” fields to be heavy since they get their masses from large VEVs 
unlike In the presence of more than — 1 fundamentals, this logic no longer applies and 

the light phase fields can constitute DM. 
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5 Dark matter phenomenology 

In what follows, we consider direct and indirect detection as well as relic abundance 
constraints on vector DM. In the U(l) and SU(2) cases, all the scattering processes are 
mediated by hi and ^ 2 . For SU(3) and larger groups, further states can contribute. 
However, we make the simplifying assumption that the Higgs mixing with those states 
is small and/or such states are heavy. In that case, it suffices to consider an exchange 
of hi and h 2 only. We note that earlier phenomenological analyses of vector dark 
matter have appeared in various contexts [I0]-[T9]. The current collider searches for 
dark matter, e.g. in the form of a monojet plus missing energy, do not set further 
useful constraints on the model, see e.g. Ref. PO] . 


5.1 U(l) dark matter 

Let us start with the U(l) case. The relevant terms in the Lagrangian are 


C D -m\A^A^ + ^ 2 ^ {—hi sin 9 + h 2 cos 9) A^A^ 

+ ^ [hf sin^ 9 — 2hih2 sin 9 cos 6* + cos^ 0) . (54) 

The couplings of hi and ^2 to the SM fields are those of the SM Higgs up to the 
suppression factors of cos6' and sin6', respectively. A phenomenological analysis of 
this model in the decoupling limit sin6' —)■ 0 can be found in ra and 

[21 HI]. Related studies have also appeared in [T21 US]- 

The DM scattering on nucleons proceeds through the t-channel exchange of hi ,2 
and leads to the following spin-independent cross section (see e.g. m), 


o-SI 

^A-N 


9^9^ m%fjf {mj^ - Av? 9 cos^ 9 

IbTT 


(55) 


where m^v is the nucleon mass and /at ~ 0.3 parametrizes the Higgs-nucleon coupling. 
One should keep in mind that there is significant uncertainty in /at and here we use 
the value somewhat smaller than the one assumed in [2]. As expected, hi and h 2 
contribute with opposite signs. Since ^ in realistic cases, the cancellation is 
not very signihcant. 



The calculation of the DM annihilation cross section is more involved due to a 
few contributing diagrams (Fig. [^. To compute the DM relic abundance, we use the 
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software package micrOMEGAs 4.1.8 |2l]. It is important to note that the leading 
s-channel annihilation amplitude due to the hi^ 2 ^exchange is proportional to the factor 


^anmh ^ COS 6 ^ ( -^-^ | , (56) 

Vs-Ki s-^hj 

where the s-parameter can be approximated by s ~ 4m^. Therefore, the amplitude is 
highly suppressed at ttia S> rnh 2 - Although other annihilation channels remain avail¬ 
able, this makes DM annihilation inefficient for heavy masses and the corresponding 
parameter space is challenged by the direct detection constraint. 

Clearly, the s-channel annihilation becomes very efficient around the resonances, 
rriA — mhj^/2 and uia — mh^/^. In this case, a very small gauge coupling is sufficient 
to obtain the right relic abundance. The hrst resonance is quite narrow due to the 
small width of the SM Higgs, whereas the second resonance is broad since many 
decay channels are available to ^ 2 . For rrih^, > the decay ^2 —)■ hihi becomes 

important. Its signihcance depends on of Eq. ([^ with increasing for 

larger (see the explicit formulae in [6]). The resonance is widened further by the 
thermal averaging over the DM momentum. 

The most important features of DM annihilation are associated with the s-channel. 
Other channels play a less signihcant role. Similar considerations apply to the indirect 
detection constraint due to the gamma ray emission in the process of DM annihilation 
(FERMI). 

The last constraint we impose is of theoretical nature. We require that the theory 
be perturbative at the TeV scale. One way to enforce it is to demand perturbative uni- 
tarity at tree level, for instance, in the process hihi —)■ hjhj. The resulting constraint 
was estimated in |2^ to be 

A, < 0(471/3) , (57) 

where A* are the scalar quartic couplings. We dehne the quartic couplings involving 
(f) by AV)juart = Xh(j,\H\'^\cj)\‘^ -|- lA<^|0|^. They can be expressed in terms of the masses, 
the gauge coupling and the mixing angle as (see e.g. m) 


Xhcf) = g sin 29 


A^ — 




Avtua 


hi 


sin^ 29 


m 


h2 


^h2 - 



(58) 


where we have used v = 2m a/9 - This implies that both Xh^ and A^ become large 
for heavy ^2 or for light dark matter. As a result, Eq. (57) imposes an important 


constraint on our model. We further require the standard perturbativity constraint 


q 


(59) 


which we hnd less signihcant for our purposes. 

Our results for U(l) DM are presented in Fig. upper row. We include the con¬ 
straints from PLANCK [23] (relic abundance), LUX [2l| (direct detection), FERMI [25] 
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Figure 2: Constraints on the gauge coupling g vs DM mass rriA for U(l) and 
SU(2). The area between the red lines is favoured by the DM relic abundance, while 
the regions above the dashed blue, dotted black and green lines are ruled out by 
direct detection, indirect detection and perturbativity of Aj, respectively. Upper row: 
U(l) dark matter with sin0 = 0.3,m/i2 = 280 GeV (left), sin6* = 0.3, ruh^ = 1 
TeV (center), sin0 = 0.2,m /12 = ^ TeV (right). Lower row: SU(2) dark matter 
with sin0 = 0.3, m /12 = ^^0 GeV (left), sin0 = 0.3, m/ij = 1 TeV (center), sin0 = 

0.2, m /12 = 3 TeV (right). 

(indirect detection) and pertnrbativity of A/. The area between the red lines is con¬ 
sistent with the thermal relic DM abundance measured by PLANCK. The LUX data 
provide the strongest constraint on the allowed parameter space. The FERMI bound 
is typically relevant for light DM, while the perturbativity bound becomes important 
for heavy ^ 2 . In each panel, the mixing angle is chosen such that, for a given m/ 12 , 
it is consistent with the LHC and EW precision data [6]. Heavier /i 2 imply smaller 
sin^, so we take sin0 = 0.3 for the left and center panels, and sin0 = 0.2 for the right 
panel. 

The two dips are associated with the resonant annihilation through hi and ^ 2 . 
The second resonance gets broader with increasing h 2 due to the increase in Xhip and 
availability of the decay ^2 —)■ hihi. The area around this resonance constitutes the 
largest parameter space consistent with all of the constraints. For m /12 = 280 GeV, 
the allowed DM mass range is about 100 GeV; for m /12 = 1 TeV, it widens to 1 TeV, 
and for m /12 = 3 TeV, it reaches more than 3 TeV. The resonance is broadened by the 
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thermal averaging over the DM momentum, so even though it appears very broad for 
mh 2 = 3 TeV, it is still consistent with perturbativity. 

The dip associated with the resonant annihilation through hi is quite narrow and 
does not open up further significant areas of parameter space. Other features of the 
PLANCK curve are local peaks corresponding to the kinematic opening of additional 
annihilation channels. For instance, the peak at ~ 80 GeV is associated with the 
W^W~ hnal state. There are further visible peaks at (m /12 + and rrih^. 

Away from the resonances, there appears to be a further allowed region at m^ > 
m/ij. Since the ^2 production is not suppressed by sin 6 *, in this case the f-channel 
annihilation AA —)■ h 2 h 2 and the quartic interactions dominate. The Planck-allowed 
strip is dangerously close to the LUX bound, so the conclusion depends strongly on 
the nucleon-Higgs coupling //v — 0.3, which suffers from substantial uncertainties. 


5.2 SU(2) dark matter 


Aside from the gauge self-interactions, the Lagrangian (54) applies to the SU(2) case 
as well, up to the summation over the 3 species, A^A^ —>■ The main change 

compared to the U(l) case is that the annihilation cross section decreases since only 
the species with the same group index can annihilate through the Higgs-like states. In 
order to keep the same relic abundance, one needs to increase the gauge coupling. Since 
the s-channel annihilation through hi ^2 often dominates, this amounts approximately 
to 

g -)■ \/lg , (60) 


whereas if the f-channel and the quartic interactions dominate, the rescaling factor 
is closer to 3^/^^. The direct detection constraint remains the same as in the U(l) 
case since particles with different group indices scatter the same way on nucleons. 
This decreases somewhat the allowed parameter space compared to the Abelian case 

(Fig. 15 . 

Non-Abelian DM features a semi-annihilation channel AA —)■ Ahi ^2 (Fig. [^. In 
some regimes, for example at large g and small sin 6 ', it can even dominate [9] (see 
also [2ni EZI)- Using the analytical results of [9], we hnd that DM semi-annihilation 
is insignificant in the relevant parameter regions. For example, at rriA > mh 2 the 
sin 0-unsuppressed and potentially important channel AA —)• Ah 2 opens up, yet it is 
dominated by A A —)■ ^ 2 / 12 . Also, around the resonances the gauge coupling is rather 
small which diminishes the relative importance of semi-annihilation. 




Figure 3: Semiannihilation of vector DM. 

No hrm conclusion can be reached as to whether the region > m /12 is allowed. 
As stated above, the uncertainties in //v play a critical role due to the proximity of 
the Planck band and the LUX bound. 
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Figure 4: Illustration of the main features of SU(3) vs SU(2) DM in the extreme 
case. Here vi/v 2 = 1-1 and sin0 = 0.3,771^,2 = 280 GeV. In the SU(3) case, Mdm 
stands for the mass of the (dominant) lighter component See Fig. for further 
details. 


5.3 SU(3) dark matter 

In the SU(3) case, DM is composed again of 3 species with two of them being degen¬ 
erate in mass and the third one being lighter (d^'). This is a result of the 

mixing between the gauge bosons corresponding to the Cartan generators of SU(3). 
Therefore, while the couplings of to the Higgs like scalars remain the same as 

in the SU(2) case, the coupling of changes. In terms of rriA = ■m^i, 2 , the relevant 
Lagrangian reads 

r 3 ( E (61) 

a=l,2 

_l_ grtiA 6 + h 2 cos 6 ) ^ ^ -f- ^ cos a — j j 

a=l,2 

+ ^ {hj sin^ 6 - 2 h,h 2 sin 0 cos 0 + hj cos^ 6 ) ( ^ + ( cos a - . 

The mass of the lighter DM component is reduced by the factor (1 — ;^tana)^/^ 
compared to that of and A^, while the gauge-scalar couplings decrease by a factor 
(cos — ^ sin a)^. Therefore, the lighter state has a smaller annihilation cross section. 

Note that sin a ~ and even a factor of two difference in the triplet VEVs leads 
to a rather small a ~ 10“^. In that case, there is no tangible difference between the 
SU(2) and the SU(3) analyses. 

To understand where differences can appear, it is instructive to consider the limit 
Vi ~ 772- Although in this case there are further relatively light states that can mediate 


18 















DM annihilation (e.g. let ns consider the simplihed example in which only the 

same states are allowed to contribnte in the SU(2) and SU(3) set-nps. The main 
featnres (Fig. are that the gange conpling mnst be larger in the SU(3) case in order 
to allow for efficient annihilation of and that the resonant dips are slightly shifted 
dne to a different freeze-ont temperatnre (see e.g. Pi). The DM density today is 
dominated by the lighter component. Since it conples to nucleons weaker than A^^'^ 
do, the direct detection constraint relaxes. Understanding further features of the 
model would require precise knowledge of the spectrum and the couplings, which we 
relegate to future work [29] . 

Finally, one should keep in mind that there exists the coupling (^3 <^ 4 , where 

and <^4 have the same parities. If <^4 were light, that is, A 4 — A 5 1 and/or vi ~ V 2 


(see Eq. (39)), the decay A^' —)■ <^94 + SM would occur. In that case, DM would consist 


of both t 


le vector and scalar components. We defer a detailed study of this scenario 


to future work 


6 Summary and conclusions 

In this paper, we have considered the possibility that the hidden sector enjoys SU(N) 
gauge symmetry and couples to the Standard Model through the Higgs portal. We 
find that when endowed with a “minimal” matter content, such hidden sectors lead 
naturally to stable vector dark matter. The underlying Lie group symmetries which 
stabilise DM are associated with complex conjugation of the group elements and dis¬ 
crete gauge transformations. 

We require complete breaking of hidden SU(N) by scalar multiplets to avoid mass¬ 
less states (barring conhnement in some cases). That can be done in a minimal fashion 
by introducing N — 1 scalar multiplets in the fundamental representation, which de¬ 
velop generic VEVs. If the scalar sector preserves CP, the above-mentioned discrete 
symmetries of the Lie group generalise to full-fledged symmetries of the model and 
lead to stable gauge fields. When sufficiently light, they constitute all of dark matter. 
In this case, DM consists of 3 components associated with an SU(2) subgroup which 
hosts the lightest gauge fields Aj^, and Hj)'. Two of them {Aj^, H^') are always 
degenerate in mass, while for iV = 2 all 3 components have the same mass. 

We have performed phenomenological analyses of U(l), SU(2) and SU(3) gauge 
held dark matter. We hnd that there are vast regions of parameter space where all 
of the relevant constraints are satished. In many of these regions, DM annihilation 
is facilitated by the broad resonances associated with the Higgs-like scalars. We also 
hnd that the SU(3) case appears very similar to that of SU(2), unless the scalar VEVs 
breaking SU(3) are close in magnitude. 
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Appendix 

A Hidden SU(3) vector—scalar couplings for U 3 = 0 

The tables below provide a list of most important gauge-scalar couplings. These are 
relevant to DM phenomenology as well as to understanding the decay channels of the 
heavier gauge helds. 


a 

b 

i 

j 

coeff. of 

4 

4 

3 

3 


5 

5 

3 

3 


6 

6 

3 

3 

8^ 1)2 

7 

7 

3 

3 


8 

8 

3 

3 

1;^2U+'"2 

6 ^ 1)2 

4 

4 

4 

4 


5 

5 

4 

4 


6 

6 

4 

4 

8 ^ v'f 

7 

7 

4 

4 


8 

8 

4 

4 

l;i2U+'"2 

6^ v'^ 

1 

4 

2 

3 

1 ~2 

2 

5 

2 

3 


3 

6 

2 

3 

1-2 a/^1+^2 

4 " Dl 

6 

8 

2 

3 

1 s2 

4a/3^ Dl 

1 

5 

2 

4 

17.2 

4 " Dl 

2 

4 

2 

4 

1-2 V^l+^2 

3 

7 

2 

4 

4^^ Dl 

7 

8 

2 

4 

1 -2^^?+’'! 
4i/3^ Dl 


a 

b 

i 

j 

coeff. of 

4 

4 

1 

1 


5 

6 

5 

6 

1 

1 

1 

1 

~gVs 

7 

7 

1 

1 


8 

8 

1 

1 

r/6 

1 

1 

2 

2 


2 

2 

2 

2 


3 

6 

3 

6 

2 

2 

2 

2 

-gVs 

7 

7 

2 

2 


8 

8 

2 

2 

~gy24 

3 

8 

2 

2 

_ ^52 

4V39 


Table 2: Non-derivative couplings 
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a 

b 

i 

coeff. of 

4 

4 

1 


5 

5 

1 

1 ~9 

6 

6 

1 

W Vl 

7 

7 

1 


8 

8 

1 

1 ~2 

1 

1 

2 


2 

2 

2 


3 

3 

2 

1 ~2 
id V2 

6 

6 

2 


7 

7 

2 


8 

8 

2 


3 

8 

2 



a 

b 

i 

coeff. of 

1 

4 

3 

l~2V2^yvf+v^ 

2 

5 

3 

id 

3 

6 

3 

1 ~2^^2^y vf+v^ 

id i;i 

6 

8 

3 

1 ~2'’2^/vf+v^ 
4a/ 3^ til 

1 

5 

4 

l~2V2y'vf+v^ 

id 

2 

4 

4 

1 ~2'"2y/vf+v^ 

3 

7 

4 

id 

7 

8 

4 

iV^d VI 


Table 3: Non-derivative couplings 
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